A CENTERLESS REPRESENTATION OF THE 
VIRASORO ALGEBRA ASSOCIATED WITH THE 
UNITARY CIRCULAR ENSEMBLE 



LUC HAINE AND DIDIER VANDERSTICHELEN 

Abstract. We consider the 2-diinensional Toda lattice tau func- 
tions Tn{t, s]r],0) deforming the probabilities Tn{rj,9) that a ran- 
domly chosen matrix from the unitary group U (n) , for the Haar 
measure, has no eigenvalues within an arc {rj, 9) of the unit circle. 
We show that these tau functions satisfy a centerless Virasoro al- 
gebra of constraints, with a boundary part in the sense of Adler, 
Shiota and van Moerbeke. As an application, we obtain a new 
derivation of a differential equation due to Tracy and Widom, sat- 
isfied by these probabilities, linking it to the Painleve VI equation. 



1. Introduction 

Consider the group U{n) of n x n unitary matrices, with the nor- 
mahzed Haar measure as a probabihty measure. The Weyl integral 
formula gives the induced density distribution on the eigenvalues of 
the matrices on the unit circle in the complex plane, and is given by 

1 " H 

-|A„(z)pr[-^; z, = e'^- and A^{z) = TT {z,-zi). 

k=l l<k<l<n 

Thus, for ri,6 e] — 7r,7r[, with rj < 6, the probability that a randomly 
chosen matrix from U{n) has no eigenvalues within an arc of circle 
(?7, 9) = {z E S^\ri < arg(2;) < 9} is given by 



^ ^ -"^ l<k<l<n 



Obviously, this probability depends only on the length 9 — rj. All of 
this is well known and we refer the reader to Mehta [1] for details. We 
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shall denote by 

R{e) = -^^\ogT4-e,e), (i.i) 

the logarithmic derivative of the probability that an arc of circle of 
length 26 contains no eigenvalues of a randomly chosen unitary matrix. 

The starting motivation for our work was to understand a differential 
equation satisfied by the function R{9) 

R{ef + 2 sin ^ cos 9 R{e)R'{e) + sin^ 6 R!{ef 

= \{\ sin^ e + sin e cos 6 R"{9) + ( cos^ 9 + sin^ 6) 

(1.2) 

obtained by Tracy and Widom in [2], from the point of view of the 
Adler-Shiota-van Moerbeke [3] approach, in terms of Virasoro con- 
straints. Introducing the 2-Toda time-dependent tau functions 

rnit,s;r,,e) = 1 / |A„(.)P J] (^^^^^'^^'^'^'^ 



(1.3) 

with = e'"^'', deforming the probabilities Tn{ri,9) = r„(0, 0; r], 6^), we 
discover that they satisfy a set of Virasoro constraints indexed by all 
integers, decoupling into a boundary-part and a time-part 

][e'''^ + e'''-^)r4t,s-,v,0)=LPr4t,s-,7^,e), k e Z, z = V^, 

with the time-dependent operators L^"^ providing a centerless repre- 
sentation of the full Virasoro algebra, see Section 2 (Theorem 2.2) for 
a precise statement and the proof of the result. 

In their study of Painleve equations satisfied (as functions of x) by 
integrals of Gessel type Eu(.n)e^ where the expectation -E;7(„) 

refers to integration with respect to the Haar measure over the whole 
of U{n), Adler and van Moerbeke P] found the s/2 subalgebra corre- 
sponding to k = —1, 0, 1, without boundary terms. The appearance of 
boundary terms and of a full centerless Virasoro algebra is to the best 
of our knowledge new. From this result, it is easy to obtain equation 
(11.21) . using the algorithmic method of [3]. Finally, similarly to a result 
by the first author and Semengue [5] on the Jacobi polynomial ensem- 
ble, we show that R{0) is the restriction to the unit circle of a function 
r{z) defined in the complex plane, so that a{z) = —i{z — l)r(z) —n^z/A 
satisfies a special case of the Okamoto- Jimbo-Miwa form of the Painleve 
VI equation. This will be explained in Section 3 of the paper. 
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2. A CENTERLESS REPRESENTATION OF THE ViRASORO ALGEBRA 

The proof of the Virasoro constraints satisfied by the integral (11. 3p is 
a non-trivial adaptation of the self-similarity argument exploited in the 
case of the Gaussian ensembles, based on the invariance of the integrals 
with respect to translations, see |6] and references therein. Here, we 
replace translations by appropriate rotations. More precisely, setting 

d/„(t,.,.) = |A„(.)r n (e^-^(*^-"+^^-'^) (2.1) 

with Za = e**^" and |A„(2;)p = ni<a</3<n — zpl'^, we have the fun- 
damental next proposition. 

Proposition 2.1. The following variational formulas hold 



dl4z^ ^ z^e<^'--^-'^) = [L^^ - L^:l) din, (2.2) 



d 
d^ 

A d/„(z„ ^ z„e-(^^+^""')) 1^^^ = + lLI) d4, (2.3) 



for all k > 0, with 



j=i dtjdtk^j dtk dtjj^k 
d 4^ . 9 



j=k+l j=l 



(n) _ _ V- "9 _ „ A _ V 7, 



^ dsjdsk-^j dsk jr[ 9sj+k 



j=i J - J - j=i 

OO r% fc — 1 o 



+ E ^'^i/^ + E^'^^/^^ + ^ ^ 1- (2-6) 

j=k+l ■' J=i ■' 

Proof. We shall only give the proof of (12. 2p . the proof of (12. Sp is similar. 
Upon setting 
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the following four relations hold, for k > 0, 



dtk 



Q = l 



(2.7) 



(9^ n 



2 ^ 2 

i+j=k ■' 
i,j>0 



--4,0 1^ 



k-1 



l<Q</3<n. 
i,j>0 



2 dsidsj 2 

i+j=k ■' 
i,j>0 



l<a<^<n 
ij>0 



(2.8) 

We split the computation into four contributions, corresponding to 
various factors in (12. ip . 



Contribution 1 : For A; > 0, we have 



d_ 
d~e 



l<a<(3<n 



le=0 



fe-l fc-1 



1<Q</3<?1 i=0 i=0 



l<a<l3<n 

2 771-1 



l<a</3<n 
i+j=k 
i,j>0 



a=l 



E. 



Using the four relations (12.71) and (12.81) . we obtain 



|-|A„(^e^(^'-^"'))|' 

I V /I 



2 771-1 



e=0 



'2\^n{z)\'E 



-E 

2 ^ 

i+j=k 
i,j>0 



dtidtj 



+ 9 E 



d"^ n — k d n — k d 



+ 



+ 



2 dsidsj 2 2 9s A: 

ij>0 



(2.9) 
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which is also trivially satisfied for /c = 0. 

Contribution 2: For k > 0, using the relations (12.7p . we have 



d_ 
d~e 



a=l 



e=0 



{{k + l)zl + {k- l)z-')E n dz. 



(2.10) 



a=l 



Contribution 3: For k > 0, using the relations (12. 7p . we have 



5_ 



n oo oo 



e=0 



-k\ 



fe— 1 n 



E 



oo n k—L n oo n 

E E 4*' - E E ' - E ^-'^ E 4-' 



fc— 1 n 



oo 71 



j=fc a=l 



- E ^'^^^ E - E E + E E 



j=l a=l j=k a=l 

oo r% fe — 1 r% oo 



E 



E 

j=k+i 



d 



nktk 



k-l 



d 



3=1 " 1 j=k+ 



CO g CO g 

E -^'^^T^r ri/ssfc + ^js,- 



1 ^Sj-fc 



9s t 



(2.11) 



Contribution 4'- For > 0, using the relations (12. 7p . we have 



d 



n 



a=l " 



e=0 



^-'[-E^^+E^.i^n^ 

a=l a=l a=l 

E-'\-^ + ^]Ef]^. (2.12) 



Adding up ([23]), (|2J0|) . (12TT]) and ( 1232|) gives (D. This concludes 
the proof of Proposition 2.1. □ 
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We are now able to state our main result. 

Theorem 2.2. (i) The tau functions^ Tn{t, s;ri,9),n > 1, defined in 
(11 -Sp . satisfy 

S,(r/, ^)r„(t, s; r], 6) = s;r],e), ke Z, (2.13) 
with Lf\ keZ, defined as m ([22D, (D, ([M]), and 

a3,(r^,^) = l(e^^^^ + e^^''|-); z = v^. (2.14) 

(ii) The operators L^k\k G satisfy the commutation relations of 
the centerless Virasoro algebra, that is 

[Lt\Lf^]={k-l)Ltl k,leZ. (2.15) 

Proof, (i) Denoting Za = e*'''", the change of variable Za ^ Zae^^^°'~^°' ^ 
in the integral (11.31) gives the following transformation on the angle 
'^a ^ fa + 2e sm{kipa), inducing a change in the limits of integration 
given by the inverse map 

(fa ^ 'fa - '^£sm{kipa) + 0{e'^), (2.16) 

for e small enough. Making the change of variable in the integral (II. 3p . 
with the corresponding change in the limits of integration, leaves it 
invariant. Thus, by differentiating the result with respect to e and 
evaluating it at e = 0, using the chain rule together with (12. 2p and 
(I2.16p . we obtain 

0= (-2sin(A;^^)^-2sin(A;77)|- + 4")-LL"i)r„(t,s;r7,e). (2.17) 
Similarly, the change of variable Z(y I y Zq/ 

^is{z^+za: ) corresponds to the 
transformation y:>Q^<fa + 26 cos{k(fa), with inverse 

which, using (12. 3p . leads to 

= ( - ^ cos(A:^)^ - ^ cos(A;r^)|- + 4"^ + L^^i) r„(t, s; r^, 9). (2.18) 

Adding and subtracting (I2.17P and (12.180 gives the constraints (I2.13p . 
with 'Bk{v,0) defined as in (ICTp . 



See the beginning of Section 3, for a justification of the terminology. 



A CENTERLESS REPRESENTATION OF THE VIRASORO ALGEBRA 7 

(ii) Consider the complex Lie algebra A given by the direct sum of 
two commuting copies of the Heisenberg algebra with bases {ha-, aj\j G 
Z} and {hf,,bj\j G Z} and defining commutation relations 

[ha, aj] = , [aj, ak] = jSj^^kK, 
[h,bj] = , [bj,bk]= j6j,^kh, (2.19) 
[ha,hb] = , [aj,bk] = , [ha,bj] = , [hb,aj]=0, 

with j, /c G Z. Let S be the space of formal power series in the variables 
ti, t2, . . . and si, S2, ■ ■ ■ , and consider the following representation of A 
in S : 



_ d_ _ . . _ A h - ■ 

o-j — 1 O'-j — jtj ) bj — ^— , o_j — jSj, 

ao = bo = fi , ha = hb = l, (2.20) 
for j > 0, and /i G C. Define the operators 

where A; G Z, a^, 6^ are as in fl2.20p with /i = n, and where the colons 
indicate normal ordering, defined by 



ajttk if j < k, 
akttj if j > k, 



and a similar definition for : bjbk :, obtained by changing the a's in 6's 
in the former. Using these notations, we can rewrite (12.41) . (12. 5p and 
(12.61) as follows 



4"' 






Tin) 










_Bin) 



^ fc-1 

- B^^l + - ^(aj - b^j){ak-j - bj^k), k > 1 



fe-i 



As shown in (see Lecture 2) the operators A^^\ G Z, provide a 
representation of the Virasoro algebra in S with central charge c = 1, 
that is 

A^"^] = {k- l)Atl + 4,-.^, (2.21) 



8 L. HAINE AND D. VANDERSTICHELEN 

for k,l G Z. Similarly, the operators -B^"^ satisfy the commutation 
relations 

, Bt^] = {k- l)Btl + 4.-.^, (2.22) 

for k,l & Z. Furthermore we have for k,l & Z 

[afc,v4{"^] = kttk+i , [&fe,5f"^] = kbk+i, 

K,i?h = , h,An = 0. (2.23) 

Let us now establish the commutation relations fl2.15p . We give the 

proof for k,l > 0, the other cases being similar. As 5 •"■'] = 0, 



i,j G Z, we have using fl^TTg]) . (EHD, and (1^^ 



Relabeling the indices in the sums, we have 



^ k+l-l 



2 

j=k+i 



^ k+l~l 



=(fc-o4i- 

This concludes the proof of Theorem 2.2. □ 
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3. The unitary circular ensemble and the Painleve VI 

EQUATION 

It is well known, see for instance [J-J, that the integral Tn{t, s]ri,6) 
defined in (11.31) can be represented as a Toeplitz determinant 

(3.1) 

with 

A nice consequence of this representation is that Tn{t, s;ri,6) is a tau 
function of a reduction of the 2-Toda lattice hierarchy, that was called 
the Toeplitz hierarchy in [4]. Therefore, as with any 2-Toda tau func- 
tion (see [8]), it satisfies the KP equation in the t = (ti,t2,---) (or 
s = (si, S2, . . .)) variables separately 

As announced in the introduction, in this section, using the method 
of |3], we establish the following result. 

Theorem 3.1. The Virasoro constraints (I2.13p . combined with the KP 
equation (13. 2p in the t variables (or the KP equation in the s variables), 
imply that the function R{9) defined in (II. ip satisfies (II. 2p . 

Proof. Remembering the definition of Lq"^ in (12. 5p . the Virasoro con- 
straint in (I2.13P for = 0, evaluated along the locus t = s = 0, gives 



d log Tn{t,s]ri, 9) 



09 



d\ogTnit,s;r],9) 



dri 

t=s=0 



(3.3) 

t=s=Q 



which is a reformulation of the fact that the gap probability r„(0, 0; rj, 9) 
only depends on the length 9 — rj. 

Define the operator D = ^ — ^ and put for a fixed n 

f{t,s;r],9) = log Tn{t,s;r], 9), 

g{r],9) = -^'D\ogTn{t,s;r],9)l^^^^. (3.4) 



Notice that for A; > 



log Tn{t,s;r],9)\t=_s=o = -r^ log Tn{t, s; -9,1 
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Clearly, from the definition of R{6) in (11.11) . we have 

Rid)=g{-d,d) = -^±\ogr4t,s;-e,e)l^^^^. 

Remembering the definition of L^"^ in fl2.4l) . the constraints in (12.1 3p 
for k = 1,2, evaluated at s = (si, S2, S3, . . . ) = (0,0,0,...), can be 
written 



S2(r/,^)/ 



s=0 



s=0 



df 


^ df 


dtj+i 


df 


+ 

s=0 


dtj+2 



(3.5) 



s=0 



+ 



s=0 



\dtj 



+ n 



s=0 



5/ 



s=0 

(3.6) 



Using (13. 3p and the definition of g{rj, 9) (13.41) . the constraint (13. 5p eval- 
uated along the locus t = s = gives 



dl 

dti 



-(e^^-e^V^,^). 
tn 



(3.7) 



t=s=0 



Consequently, along the locus r] = —9, we have 



dl 
dti 



t=s=0 
r}=-e 



sin(^)i?(^). 

n 



We then proceed by induction. We call 

Quf 



dtj,dtj,...dtj^ 

a t derivative of weighted degree \j\ = ji + J2 + ■ ■ ■ + jn- Then, for 
> 1, we compute the system formed by 

all t-derivatives of weighted degree k of (13. 5p . 
all t-derivatives of weighted degree A; — 1 of (13.61) . 

evaluated at t = s = 0. For instance, for = 1, (13. 8p reduces to 



(3.8) 
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After substitution of (13. 7p . this system of equations can be solved for 



^ ^^^^Q can be expressed in terms of 6, R{6) and R!{6). 

ri~=-e 

For general > 1, suppose all the t-derivatives of / of weighted 
degree fc, evaluated at t = s = 0, have been expressed in terms of rj, 6 
and g{ri, 9), ... ,'D^~^g{ri, 9), whenever n ^ 1, . . . ,k — l. Then (13.81) is a 
system of linear equations where the unknowns are all the t-derivatives 
of / of weighted degree k+1, evaluated at t = s = 0, and the coefficients 
can be expressed in terms of t], 9 and g{ri, 9), ... ,'D''~^g{r], 9). This is 
a system of p{k) + p{k — 1) linear equations in p{k + 1) unknowns, 
where p{k) is the number of partitions of the natural number k. As 
p{k + 1) < p{k) + p{k — 1), this system can be solved and all the t- 
derivatives of / of weighted degree + 1, evaluated at t = s = can be 
expressed in terms of rj, 9, and g{ri, 9), . . . , 'D^g{ri, 9), whenever n ^ k. 
Consequently, on the locus r] = —9, the t-derivatives of / of weighted 
degree k + 1, evaluated at t = s = and on the locus r] = —9, can be 
expressed in terms of 9, R{9), R'{9), R^''\9). 

Since the KP equation (13.21) contains t-derivatives of / of weighted 
degree less or equal to 4, by performing the above scheme up to A; = 3, 
we can express all these derivatives, evaluated at t = s = and r] = —9, 
in terms of 9, R{9) and its first three derivatives, whenever n > 4. This 
gives us a third order differential equation for R{9): 



= 4R{9y - 2(^2 + (1 - n^) cos29)R'{9) + 8 sin 29 R{9)R' (9) 

-2sm29R"{9) + sin^ 9 {l2R! {9f - R"'{9)). 

Multiplying the left-hand and the right-hand side of this equation with 
i sin 9 h cos 9 R\9) + sin 9 R"{9)) , we obtain 





(3.9) 



with 



W{9) = R{9f + 2 sin ^ cos 9 R{9)R'{9) + sin^ 9 R'{9f 
--(\ sin^ 9 + sin 9 cos 9 R"{9) + ( cos^ 9 + 7? sin^ 9) R'{9)) 
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Equation (13.91) implies that W{9) = 0, which is the equation (11. 2p . 
obtained by Tracy and Widom in This concludes the proof of 
Theorem 3.1. □ 

Remark 3.2. In the above proof, we had to assume that n > 4, where 
n is the size of the random unitary matrices. For n = 1,2, 3, the func- 
tion R{6) also satisfies (11. 2p . as can be shown by direct computation, 
using the representation (13. ip of the probability t„(?7, 9) as a Toeplitz 
determinant. It would be interesting to relate the proof with the origi- 
nal derivation in [2] • For the Gaussian ensembles, the relation between 
the two methods has been studied in [9] . 

Finally, similarly to the case of the Jacobi polynomial ensemble (see 
[5]), we observe that R{9) in (11.11) is linked to the Painleve VI equation. 
Precisely, we show that it is the restriction to the unit circle of a solution 
of (a special case of) the Painleve VI equation, defined for z E C 

Corollary 3.3. Put R{9) = r(e~2i^). Then, the function 

a{z) = -i{z - l)r{z) - —z 

satisfies the Okamoto-Jimho-Miwa form of the Painleve VI equation 

[z{z - l)a"f + Az{z - l){a'f + Aa'a^ + 4(1 - 2z)a{a'f 

- ci{a'f + [2(1 - 2z)ci - C2](t' + Ac^a - Cg = 0, (3.10) 

with 

o 4 6 4 

2 on n n 

ci=n, C2 = — , C3 = — , C4 = - — . (3.11) 
8 lo lo 

Proof. From (II. 2p . by a straightforward computation, putting R{9) = 
r(e~^*^), we obtain that r{z) satisfies 

[z{z - lyf + Az\z - l)r'r" - Aiz{z - l^ir'f - At{z^ - l)r(r')' 

+ [4z^ - n\z - lf]{r'f - Airh' = 0. (3.12) 

Substituting in (I3.12p 

, , aiz) + xz 
r[z) = I — 

for some constant x, and annihilating the coefficient of cr^, one finds 
that X = n^/4. With this choice of x, the new function a{z) satisfies 
the Painleve VI equation (I3.10p if we pick ci, C2, C3 and C4 as in (13. lip , 
which establishes Corollary 3.3. □ 
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4. Discussion of the results and some further directions 

Our starting motivation was to understand a differential equation 
(11 ■2p due to Tracy and Widom [2J, satisfied by the logarithmic deriva- 
tive of the gap probability that an arc of circle of length 29 contains no 
eigenvalues of a randomly chosen unitary n x n matrix, from the point 
of view of the algebraic approach initiated by Adler, Shiota and van 
Moerbeke [3]. The main surprise is that the 2-dimensional Toda tau 
functions (ll.3p deforming these gap probabilities, satisfy a centerless 
full Virasoro algebra of constraints. The result stands in contrast with 
the corresponding integrals for the Gaussian or the orthogonal polyno- 
mial ensembles, which roughly satisfy only "half of a Virasoro type 
algebra of constraints, see [3], [5], [6] and 

As mentioned at the beginning of Section 3, the integrals (II. 3p can 
be expressed as Toeplitz determinants, see (13. ip . As such, they are very 
special instances of tau functions for the so-called Toeplitz lattices [1], 
that is 



Tn{t, s) — det(/ifc_/(t, '5))q<^ 

(4.1) 



where 



IJk(t,s)= [ z'' e^T^^^'^''+'^'~'^w(z) fceZ, (4.2) 

and w{z) is some (complex- valued) weight function defined on the unit 
circle S^, such that the trigonometric moments 

f dz 

/ifc = /ifc(o,o) = / z''w{z) ——, kez, 

Jsi 2mz 

satisfy det [fj^k-ij^^f^ ;<n-i 7^ 0, > 1. In the special case (13. ip that we 
consider in this paper, w{z) = X(r;,e)=(-2) is the characteristic function of 
the complement of the arc of circle [i], 6) = {z & S^\ri < aig{z) < 6}. 

As it immediately follows from (14. 2 p , at the level of the trigonometric 
moments, the Toeplitz hierarchy is given by the simple equations 

_ (9/ifc _ _dfik _ w ■ ^ T 

Obviously [Ti,Tj] = 0,^i,j G Z, if we define To/i^ = fj,k- Following an 
idea introduced in |S] in the context of the 1-dimensional Toda lattices, 
we define the following vector fields on the trigonometric moments 



Vjfj,k = {k + j)fJ,k+j, Vj G Z. 



(4.3) 
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These vector fields trivially satisfy the commutation relations 

[V„V,] = ij-t)V.+, (4.4) 

[V„T,]=jT,+„ Vz,jGZ, (4.5) 

from which it follows that 

m,T,],T^] = j[T,+„T,-] = 0, J G Z. (4.6) 

Equations (14. 4p . (14. 5 p and (14. 6 p mean that the vector fields Vj,j G 
Z, form a Virasoro algebra of master symmetries, in the sense of 
Fuchssteiner pIJ], for the Toeplitz hierarchy. 

The tau functions (14. ip admit the following expansion 

'Tn{t,s) = ^ Pio,-,jn-l'S'i„_i_(n-l),...,jo(^)'S'j„_i-(n-l),...,jo('5)) 

US*0<'"<«n-l 

o<io<---<in-i 

where 

Pio,...,in-i = det (/ii,_j-,(0,0))(,^^^,^^_^, (4.7) 

are the so-called Pliicker coordinates, and Si-^^„_^if,{t) denote the Schur 
polynomials 

Sii,...,ik(t) = det (S'j^+s_r(t))^^^ 

with ^^(t) the so-called elementary Schur polynomials defined by the 
generating function 

oo 

exp ( ^kX^) = ^nitiM. ■ ■ Oa;"- 
k=l neZ 

In a forthcoming publication [11], we shall establish the next result: 
Theorem 4.1. For all k & Z, we have 

Vkipio,---,in-l)Si„_^^(n-l),...,io(t)Sj„_^-(n-l),...,jo{s), 
0<^0<■■■<i„-^ JO,...Jn-i^ 

o<io<---<in-i 

with L^^\ A; G Z, defined as in (12.41) . (12. 5p . (12.61) . and Vk(pio,...,i„_i ] the 

Lie derivative of the Pliicker coordinates (14. 7p m i/ie direction of the 
master symmetries Vk of the Toeplitz hierarchy, as defined in (14. 3p . 

Thus the operators L^"^ , A; G Z, precisely describe the master sym- 
metries of the Toeplitz hierarchy on the tau functions of this hierarchy. 
Since master symmetries are usually connected with a bi-hamiltonian 
structure in the sense of Magri ^12| (see [13] and [H] for an overview). 
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it suggests to investigate the existence of a bi-hamiltonian structure for 
the Toephtz hierarchy, which seems to be an open problem. 
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